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_ 3 . ABSTRACT 

An approximate analytical queuing model for the intraship transfer 
of material phase of the replenishment at sea problem is formulated. 

The physical model has two parallel sets of three service facilities 
in series that move material to a final stage of three service 
facilities in parallel. Exponential service times and finite queue 
space are assumed for each facility, except that the initial facilities 
always have infinite queues. This implies that the output from an 
individual facility is approximately Poisson. The mean output rate for 
the queuing system is obtained. Several parametric studies on critical 
parameters are performed, with the results presented in graphical format. 
The model is easily enlarged to allow a general number of facilities in 
series in any of the parallel paths. Generalization to more than two 
paths in parallel for the initial stages may be possible. 
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I. INTRODUCTION 



Underway replenishment (UNREP) of the ships of the United States 
Navy at sea is a mandatory requirement in maintaining a continuous 
posture of readiness to meet the present and future global commitments 
and responsibilities of this country. However, replenishment of ships 
at sea should detract as little as possible from the primary missions 
of a fleet. This implies that the goal of the UNREP operation is 
complete resupply of combatant ships in the minimum possible time. 

There are several phases of the replenishment operation , ^ possibly 
not all distinct. One list of phases might be: Long range planning 

required to schedule a replenishment both in terms of time and place; 
the actual movement of ships toward the prescheduled rendezvous; when 
several supply and combatant ships are involved in a replenishment 
operation, the scheduling and movement of all the combatant ships 
through the replenishment cycle; and the actual transfer of material 
from the hold of the supply ship to designated storerooms on the 
combatant. It is obvious that all phases require careful attention 
if the operation is to be efficient. Consequently, all phases are 
fruitful study areas in a continuing effort to improve UNREP efficiency. 

The study reported herein dealt with the actual transfer operation. 
Perhaps it is the most important phase in that the ships involved are 
most vulnerable, their capability to maneuver reduced and their 



Throughout the paper the terms UNREP and replenishment operation 
will imply the entire operation; planning, scheduling, shipmovement , 
etc. , while the term transfer operation denotes the movement of material 
from hold of the supply ship to storeroom of the combatant. 
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watertight boundaries open to allow passage of material. Unfortunately 
in the past it appears that this phase has received less attention than 
it warrants both from the standpoint of detailed study and the introduc- 
tion of improved material handling equipment and techniques. 

Recent efforts have been made that should considerably improve the 
overall speed and efficiency of the transfer operation. Improvements 
include the introduction of new supply ships with the capability to 
employ the unit load system of material storage and transfer, more and 
better mechanized material handling equipment on both supply and 
combatant ships, and a realization of the problem by ship designers 
which has resulted in improved match of replenishment stations between 
supply and combatant ships. The unit load concept provides for mate- 
rial to be stored on the supply ship in the configuration for future 
transfer. Full implementation of the concept will maintain the unit 
load intact through ultimate storage in the storeroom or magazine of 
the combatant receiver. Although this concept may never be completely 
feasible, especially in the case of smaller combatant ships, it should 
provide a vast improvement to perhaps the most important and demanding 
of all UNREP operations, the rearming of aircraft carriers. 

The transfer operation can be represented as a queuing problem. 

It is a system of parallel queues in scries with movement between the 
parallel paths. Attempts at reaching an analytic solution are further 
thwarted by the finite waiting room in the queues involved. Formal 
efforts to analyze the problem have been few in number. The powerful 
computer simulation tool has been the general mode of analysis. An 
early simulation was conducted by Stanford Research Institute [fiT] . A 
detailed computer simulation study has also been done by Alfred M. 

Fcilcr of Project TRANSIM, University of California at Los Angeles, 
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Department of Engineering. Computer simulation studies are currently 
being done at Stanford Research Institute and Center for Naval 
Analyses . 



II. BACKGROUND 



A. THE TRANSFER OPERATION 

Although some transfer of material between ships at sea is 
presently accomplished by helicopter, the vast majority of material 
moves between ships via attached horizontal rigs. A typical move- 
ment sequence for material transferred at sea would commence with a 
vertical lift from the hold of the supply ship to the main or transfer 
deck by means of an elevator or a boom and windlass. Once on the 
transfer deck the material will be moved to the transfer station by 
fork lift truck or hand pallet truck. The intraship transfer of 
material is accomplished by one of several different cable and pulley 
arrangements. The type employed on a given transfer will depend upon 
the types of ships involved and the amount and type of material to be 
transferred. Material landed on the transfer station of the combatant 
ship is cleared from the area by fork lift or pallet truck. Eventually 
the transferred material will be transported vertically down to its 
designated storeroom or magazine by elevator. 

The type of movement just discussed employed the unit load concept 
required by the model developed herein. As noted earlier, many small 
combatant ships must break up the unit loads into smaller units for 
movement beyond the transfer station. Many times hand movement is 
required because of space limitations in the ship's internal passages. 
The approximate model developed by this study could still be used for 
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this type of movement using a mean service time for the unit load 
transferred by pieces. 

Queues of material and subsequent blockages of material can and 
do form at several points along this path. On the supply ship, a 
queue of material will probably exist on the transfer deck near the 
point where vertical movement from the hold is completed. A queue 
of unit loads will also exist at the transfer station as material 
awaits movement to the combatant. The transfer station on the 
receiving ship is a critical blocking point as transfer of a load 
from the supply ship to the combatant may not begin until the landing 
area has been cleared of the previous load. Generally, material 
cleared from the landing area on the combatant ship is temporarily 
stored in what is referred to as a marshalling area. Material in 
the marshalling area is distributed to a strikedown route and then 
to eventual storage in the assigned storeroom. In many cases, 
additional queues of material will form along this strikedown route, 
especially at the top and bottom of vertical movement leg(s). 

Much of the material transferred at sea presently is stored as an 
individual component or packaged with similar items in such a manner 
that hand movement is possible. These individual components are 
combined into loads that can be conveniently handled by the transfer 
rig between the two ships. This load makeup is accomplished either 
in the hold or on the main deck of the supply ship. The process is 
essentially reversed on the receiving ship with the breakup of loads 
commenced early in the storage routine, especially on smaller ships, 
where it is generally cleared from the landing area piece by piece. 
Introduction of more and better handling equipment on all ships and 
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movement of material via unit loads through all or most of the cycle 
will vastly improve the efficiency of the operation. This step forward, 
although far from universal reality, is in progress. 

Careful consideration indicates that the three distinct entities 
having a major effect on the movement of material through the system 
are the service or handling rates through all phases, the storage 
capacity of the various queues that inherently form and the origination 
and final distribution of the transferred material. The importance of 
the handling rates is obvious to the extent that discussion here is un- 
necessary. However, the importance of the queuing (temporary storage) 
space and the distribution of material requires some explanation. Any 
amount of free open space on any ship is always scarce, this includes 
space for temporary storage of transferred material. Perhaps the single 
exception is the hanger deck on an aircraft carrier, (aircraft carrier 
transfer stations are located off the hanger deck), however, even here 
space may be at a premium as it is normally taken up by aircraft in 
various states of overhaul. Therefore, the transfer process is often 
slowed because of blockages caused by the inherent space limitations 
on board ship. 

Limited space on ships is also partially responsible for transfer 
slowdowns that can be attributed to excessive distribution to a single 
material handling facility, e.g. , a vertical strikedown route on the 
combatant ship. Space available on a ship restricts both the size of 
storerooms and vertical and horizontal passages leading to the store- 
rooms. In order to provide for convenient and rapid retrieval of a 
given item, it must be stored in a predesignated place. Further, 
many items require special storage facilities that are available in 
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only a single location. These special facilities £an range from 
refrigerated storage for fresh and frozen provisions to storage trays 
and racks for most ordnance. Even if the above difficulties could be 
overcome, space limitations again would be limiting in that storerooms 
are generally kept full so that the only storage space available is 
the designated one. Therefore, it is often impossible to change 
material destination to accommodate transfer efficiency. A possible 
solution to blockage on a route would be an alternate route for some 
of the material. Unfortunately, alternate routes may not always exist 

B. PARALLEL/SERIES QUEUING MODELS 

Solutions for a single stage of several parallel queues, with 
various types of input and service times, have been used for some time 
Knowledge of single server queues in series is more limited and is 
generally restricted to Poisson input to the first stage and exponen- 
tial service times in all stages with infinite queues allowed. 
Information in the literature of a combination of queues in series 
and parallel is virtually nonexistent. 

Saaty [6j gives a survey in the area of series queues. The 
bibliography for this work is nearly complete and he has recently 
addended the bibliography to include recent contributions in the 



fie Id 7 . Burke 1 reported the fundamental result for queues in 



series with Poisson input to the first stage, exponential service 
time in all stages with infinite queues allowed, showing that the 
output from each stage and therefore the input to the next stage 



times did initial work in the area of allowing only finite queues to 



is also Poisson. Hunt using the assumption of exponential service 




form before each series facility. 




extended 
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Hunt's work under the exponential service time assumption and also 
solved the problem for Erlang service times for each stage. 

Hillier and Boling concurrently reported an approximate procedure 
for determining the mean output rate for a system of exponential 
service time facilities in series with finite queues allowed before 
each facility except the initial stage where an infinite queue was 
assumed to exist. Hillier and Boling compared results obtained with 
the approximate procedure with the limited exact results obtained. 

The comparison showed very little error in all cases with the magnitude 
of the error decreasing with increased allowable queue size. 

The Hillier and Boling approximate procedure has the distinct 
advantage of being computationally feasible, even for very large 
problems. This approximate procedure formed the basis for the study 
reported here and it will be discussed and outlined in some detail in 
the next section. 

C . APPROXIMATE PROCEDURE 

The important features of the Hillier and Boling procedure, 

for approximating efficiently the mean output rate (R) for a series 
of single server exponential service time facilities, may be outlined 
as discussed below. Assume that there are N facilities in series each 
with an exponential service rate ofy^, i = 1,2, *** , N, and waiting 
room M^, i = 2,3, *** ,N, where the M^ ' s include the item being served 
and that an infinite queue always exists before the initial facility. 

The basis for the procedure comes from the results due to Burke [lj 
which indicated that when the M^'s are equal to infinity then each 
service facility would have a Poisson input. Reich [ 5 ] showed further 
that in this case the queue sizes would be independent. This implies 
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that each facility could be analyzed separately using the standard 

M/M/1 model. The approximate procedure with finite queues allowed 

is then based on the assumption that the output from each stage is 

approximately Poisson. Each facility is analyzed individually by 

using the single-server queuing model M/M/1 with finite queue allowed. 

Note that R is merely the mean output rate for each of the facilities. 

The procedure then assumes that there exists an effective mean 

arrival rate /\(j) (when there isn't blocking in the (j-l) St facility) 
ef f 

to the j facility. It is also assumed that the queue lengths before 
each facility are independent so that the unconditional probability 
(1-P^j - ^) that the (j-l) st facility is busy is equal to the conditional 
probability that it is busy, given that the queue for the j facility 
is not full. Thus it follows that, 



(1) ^ff = •0/j- 1 a ' P o J ’ 1) ) • 






o 



, if n = Mj+1 



where n is the number of customers in the queue for the j 1 -* 1 facility 

s t 

(including any customer being held at the (j-1) facility because of 
blocking) and j = 2,3, , N. 

The effective service rate for the j 1 "* 1 facility, will be 

somewhat less than j = 2 , 3 , ,N-1 , because time spent in holding 

a customer after service due to blockage in the (j+l) st queue must be 
included. Now the mean output rate of a single- server queuing system 
is the product of the mean service rate and the probability that the 
server is busy. Therefore, it follows that, 



(2) ^ ( J f \ = R/(1-P^ j) ) , for j = 2,3, *•* ,N. 

However, for the final server the customer is always released immedi- 
ately, i.e.,^^^ = . Therefore, 
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( 3 ) 




Thus R may be determined if is known. 



o 



For a single server queuing system with finite queue, it can be 
shown that 



(4) 




where for the case of the approximate system, 





It is also clear that since the queue before the first facility is 



can be determined from equation (1) . It then can be used along with R* 



calculated from (3) . This calculated value of R is compared to R* and 
a new trial value assumed and the process repeated until the true value 
of R is obtained when R* - R = 0. 

Hillier and Boling have simplified the above procedure for solution 
by numerical methods and the equations used in this study are presented 
below: 

Combining equations (1) and (2) yields, 



infinite, P (1) = 0 or (1-P (1) ) = 1. 
’ o o 



It is now possible to determine the true value of R by first 
assuming a trial value of R, call it R*. Since (1 -Pq^) is known, 




in equations (2) and (4) to determine (1-P^)). Then l-P^-^), 

••• , (1-P^^) can be obtained in a similar manner so that R can be 

’ o 





(j_1) )(l-P (:i) )/R , j = 2,3, ••• , N-l 



o 
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then , 



(5) 



1-P 



O) _ 



(j-D 

^j-l (1 - P o > 



y ; , j = 2,3, ••• , N-l, 



so that y can now be found by combining (4) and (5), as the positive 
real root of the equation, 



M. + l 

(6) (1-y.J ) 

y j * 



R 



y. = 0. 



M.+2 ,, _(j-l) 

ly j J Xj-i (, ' p o > 



Hillier and Boling show that this root exists and is unique for values 
of R that are sufficiently close to R, the true mean output rate. 

Let , 

M. + l 

v . (1-y . J ) 

f (y ) = J J » j " 2,3, *" ,N-1 

j j , M.+2 

i - y J 



and 



g.(R,y ) = f.(y.)-c y , j = 2,3, ••• , N- 1 
J j J J J 



where R is a trial value, 



c . = 



R 



j ~ M j-i x I-i 



* * * 

and = fj(y^), y j corresponding to the unique root of g(R,y^) 



= 0, 



Finally X yields the calculated value of R where 



< ■ f N <Aszl Cl>- 

s* N 
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III. FORMULATION OF THE PROBLEM 



For the replenishment at sea transfer problem, as with many other 
physical problems that can be formulated in a queuing' theory context, 
one of the interesting questions that can be asked is, "How long will 
the operation take?" This question was answered in the formulation 
for the UNREP transfer problem by determining the mean output rate (R) 
for the system. 

As previously noted, each at sea transfer of material will be 
somewhat different from a queuing standpoint because of the ships 
involved, volume and type of material transferred, material origin 
and destination, types of material handling devices employed, etc. 
Thus, the number and size of queues formed during any transfer may 
differ considerably. 

A completely general model to represent the transfer at sea opera- 
tion would require up to five parallel paths from the supply ship to 
the combatant ship. On each of the parallel paths there are several 
material handling facilities with associated queues and waiting space 
on both the sending and receiving ships. Additionally, the general 
model would require movement of material from one parallel path to 
another at least once on each ship. 

The general model described above may be simplified considerably 
and still represent accurately many transfer operations. The model 
developed for this study utilized two parallel paths from the hold 
of the supply ship to a marshalling area on the combatant ship. From 
each of the two marshalling areas the material is moved to any one of 
three strikedown routes. Figure 1 is a pictorial representation of 
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this model. The inclusion of only two parallel paths was done for 
computational efficiency, noting also that the vast majority of at sea 
transfers are conducted with but two attached rigs. A generalization 
of the model for more than two paths is discussed in section VI. 
Restriction of movement between parallel paths on the supply ship is 
the other major simplification made for this study. It was felt that 
such an assumption was not overly restrictive and accurately represents 
many actual transfers. 




Figure 1 



If individual material handling facilities in parallel are referred 
to as a stage, then the physical location of the first stage queue of 
the model can represent the hold of the supply ship. The second stage 
queues are formed by material on the main deck of the supply ship 
awaiting intraship transfer. The marshalling areas, which will gener- 
ally be adjacent to the transfer stations on the receiving ship, form 
the waiting space for the third stage queues. The fourth stage queues 
may be located physically at any one of several places on the individual 
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strikedown routes to the designated storerooms, e.g., on the transfer 
deck near an elevator or hatch, or at the bottom of a vertical movement 
leg. 

As stated previously, service of material in the two third stage 
queues consists of movement to any one of the three strikedown routes. 

In the model, this movement is governed by two probability distribu- 
tions p^ and q^, i = 1,2,3, for the upper and lower routes respectively. 
The distributions are determined as that percent of the material moving 
across one of the initial routes that is scheduled for eventual storage 
in a storeroom on the i*"^ strikedown route. Some of the factors that 
effect this determination are the amount and type of material trans- 
ferred, the storage of material on the supply ship, the location of 
the storage space available on the combatant ship, and the availability 
of alternate routes to final storage. It is felt that the distribution 
approach as a means of representing the actual movement during transfer 
is reasonably flexible to cover accurately many transfer operations. A 
primary consideration in this regard is the fact that material is 
selected for movement in the most convenient manner, i.e., a load is 
selected for movement because it is available to the transfer facility 
or because it can be moved ahead of other unit loads to a nonblocked 
facility upstream. 

Within the model it is assumed that all material is moved as a 
unit load. Solution for smaller units would overly complicate the 
problem, however, the unit load assumption is enhanced by the present 
trend toward actual movement of material in this manner. This assump- 
tion had a further and probably a greater impact on obtaining an 
eventual solution because of the effect it had on the transfer rates. 
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Since all material moves in this manner it is reasonable to assume that 



a given facility can handle a unit load of 500 lb. bombs or 1000 lb. 
bombs or even fresh provisions, for example, at approximately the same 
rate. Thus the unit load concept helps justify the requirement that 
each facility must move all types of material at a fixed rate during 
each transfer investigation. There is of course some allowance for 
variation, e.g., if two types of material had different service rates 
for a given stage and each moved on different parallel paths, then the 
service rates for the two paths could be adjusted to represent more 
accurately the material being handled. 

Notation within the model is such that the first subscript (i) 
represents the i 1 "* 1 stage in series and the second subscript (j) repre- 



rate for the second parallel facility in the first stage and is the 

number of unit loads that can queue before the first parallel facility 
in the third stage. 

IV. PARALLEL SERIES SOLUTION 

The flow of material through the physical model of figure l may be 
represented as in figure 2. If R is the mean output rate, it follows 



sents the j parallel element in a stage, e.g.. 




the service 



R 



R 



2 



l 




J 



R 



Figure 2 
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from previous discussion of the p.. and q.. distributions and 
figure 2 that , 



R - R. + R 
1 2 



R Rp + Rjj + Rj-jj 



R I = Pl R l + ^2 R 2 



R = p R + q R 
II 21 22 



R III P 3 R 1 + q 3 R 2 



where R^, R^, R^. , R^, and R ate the mean output rates of their 
respective routes. For any UNREP that can be represented by this model, 
the value of the mean output rates will be determined by the various 
service rates, queue waiting space, and the p^ and q^ distributions. 

If the system were symetric, that is if the service rates and queue 
lengths for the facilities in each stage were the same and the p^ and 
q^ distributions moved material symetrically to the final stage, then 
R would equal R^. A solution to this problem could be easily reached 
by the Hillier and Boling procedure by considering only half the system 
and solving for R/2. However, these conditions will not often hold and 
in general, R will not equal R^. 

When R^ and R^ are not equal, the Hillier and Boling approximate 
procedure cannot be applied directly to the system to obtain a solution 
for R. The problem is that not only must a trial value for R be 
assumed, but also values for R^ and R^ must be assumed such that 
R ! + R 2 = w ^ ere ' s are trial values. Now there are infi- 

nitely many combinations of R = R* + R* - R ,f . The problem thus 
reduces to finding the values of Rp and R 2 that will yield the true 
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mean output rate. The true values of the R's will be complicated 



functions of // , M , p , and q . It will not, however, be necessarv 

^ij ij i i 

to determine this function, as the solution may be reached by simply 

maximizing the sum + R 2 for a given value of R". 

Before proceeding to show how the true values of R^, R^ and R may 

be obtained, the Hillier and Boling equations should be adapted to the 

problem. Assume that for a given problem the true values of R^, R^ , 

and R are known. The solution for the output across the upper route 

may be found by applying the approximate procedure previously outlined. 
★ 

First may be found as the unique positive real root of 

ic 

g (R.y ) = f (y )-c y =0 
6 2i 1 21 21 21 2121 



where 



> -<A(i x u • (x n = 1K Then ’ 






,, * M 9 1+k 

y (l-y 21 ) 

V - r / x 21 21 

X 21 ■ f 2 1 ( V 2 1 \ *M +1 

21 

y^^ and then X^j may be found in a similar manner. Also, ma y be 

found using the same technique using R^ in place of R*. Then since 
the fourth stage is the final stage X,. i * 1,2,3, can be determined 
directly as 



X* = f 
4 i 4 i 



( ^l X 31 P i + ^32 X 32 q i), i = 1,2,3 

^4 i 



so that finally R = \ x + + * . 



The solution procedure for the interesting case of unequal parallel 
outputs can now be outlined. As before, R is assumed and them R^ and 
R 2 are chosen such that R^ + R 2 = R*. The problem is then solved as 

just outlined to obtain a computed value of the mean output rate, R. 
Using the same value of R , R{ and R 2 are varied (maintaining the 
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relationship = R ') until the maximum computed value of R is 

found.) As in the original Hillier and Boling procedure, the computed 



value of R is compared to R“ and if R = R“ , the true mean output rate 
has been found. If R is not equal to R , then a new value of R' is 
assumed and the procedure repeated until the true mean output rate is 
obtained. This procedure converges rapidly if two techniques are 
employed. When selecting a new value for R“ , select R x ( new ) = (R % +R) /2 . 
Also when commencing each iteration R£ (new) '^2(new) = Rl /R 2 . 

Justification for this procedure follows from physical flow laws. 

It is easy to see that the procedure is correct at the boundaries. For 
example, assume that the system was such that one of the j ' s = 0 in 
the lower path. Then maximum R will be obtained when R = R^ and R£ = 0. 

A computer program to solve the problem of figure 1 is presented on 
page 34. The program was employed to obtain the results presented in 
section V. The wide range of input data was handled with no difficulty. 
Average computation time for the IBM 360 was 2 seconds. 



V . PRESENTAT ION O F RESULTS 

As noted previously, the three types of parameters within the 
approximate transfer model that have an effect on the mean output rate 
are the service rates of the material handling facilities (y^Y^j’s), the 
available waiting room in each queue (M^j's), and the distribution of 
material between the series facilities as governed by the p^ and q^ 
distributions. The model can be employed to study UNREP transfer 
operations in at least two ways. First of all, general parametric 
studies can be made to determine the effect of varying the parameters 
ir all possible combi nations . Also, more specific studies may be made 
employing known handling rates, queue space and expected material 
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distribution for specific ship types, to determine critical areas that 
can be improved to increase overall transfer efficiency. A second use 
for the model would be as an aid to the commander, perhaps even on as 
low a level as commanding officer of a supply ship, to help him in the 
planning for replenishment operations. All that would be required is 
a knowledge of material handling characteristics of the ships involved, 
the amount of material to be transferred and a small computer which is 
generally carried on board many supply ships today. 

The primary advantages of the approximate procedure over known 
computer simulation models are the simplicity of the analytical model, 
the inherent computing efficiency and the ease with which inputs into 
a computer program may be changed. Possible objections to the model 
over computer simulation models would include the limitation that 
different types of material must be transferred through a facility at 
the same rate and that unit loads of each type of material would occupy 
the same volume or area of queue space. The exponential service time 
assumption with its memoryless property, i.e., the probability of a 
unit completing service in a time interval At is constant, no matter 
how much time has elapsed since the last completion of service, also 
may be partially objectionable. A final objection may be raised over 
the final output of the model which is a steady state result, i.e., it 
is assumed that the system has been in operation for a long time and 
is in equilibrium. This objection is at least partially overcome by 
the generally large volume of material transferred during the UNREP's 
that are of interest and the comparatively small queue sizes allowed, 
thus assuring that equilibrium conditions will be nearly attained 
during an operation. 
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Actual UNREP transfer data was not available in the format required 
for entry in the model. Therefore, numerical results were necessarily 
confined to several parametric studies. Graphs of the study results, 
along with the input data for the constant parameters, is presented in 

f 

Appendix A. 

Parametric studies 1, 2 and 3 all vary the parameter p^, distribu- 
tion of material from route to route R^. , over the range 0.333 to 1.0. 
The parametric variation in the three studies involves the distribu- 
tion as follows: 

Study one (page 29) - - - = 1/3, i = 1,2,3 

Study two (page 30) - - - q^ = 0.2, q^ = q^ = 0.4 

Study three (page 31) - - q^ = 0.4, q^ = q^ = 0.3 

All three studies indicate the marked decrease in overall transfer 

rate as congestion at the upper strikedown route (facility (4,1)) 

increases. As predictable, the decrease is most rapid for study 3 

where the flow to facility (4,1) and, therefore, the congestion at the 

facility is greatest. An additional effect, evident in all three studies, 

is the difference between R, and R that can be attributed to the slight 

1 2 

reduction in handling rates of facilities (2,2) and (3,2) and the 
smaller queue allowed at facility (3,2). 

Parametric study 4 gives a good indication of the efficiency that 
can be gained by increasing the handling rates of the facilities in a 
single stage. In this case, the service rates of facilities (2,1) and 
(2,2) were varied together over the range 0.5 to 1.0 loads/min. with 
all other service rates fixed at 0.5 loads/min. The graph of the study 
(page 32) indicates that small increases above the system average will 
have a significant effect. Study 5 varied the queue space in facilities 
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.(3,1) and (3,2) over the range 5 to 10 loads with the remaining queue 
spaces held fixed at 5 unit loads. Once again a graph of the results 
(page 33) indicates that small changes above the norm are worthwhile. 
Comparison of studies 4 and 5 shows that increased service rates have 
a greater effect on efficiency than increases in allowable queue size 
This is understandable and also encouraging, as generally, space on 
ships will not be available to increase temporary storage to promote 
increased transfer efficiency. 



VI. EXTENSION OF THE MODEL 



There are several ways that the specific model developed here can 
be extended. Unfortunately, the most interesting extension becomes 
computationally difficult. 

The most obvious extension follows directly from the Hillier and 
Boling j_2^] series development where the number of stages in series is 
N. The parallel series model developed here may be similarly extended. 
Referring to figure 2, this would mean that there could be N^ and N^ 
separate service facilities in the two initial parallel routes and , 

N , and N^^. separate service facilities in the three final routes, 
where N^ , , N^ , N^ ^ , and N^.^ are only restricted to be finite. 

Solution of this generalization should present no special difficulty. 

The interesting extension fron an UNREP point of view would seem to 
be that which generalizes the number of parallel paths allowed. Such a 
generalization is represented in figure 3. 

In solving such a system, it is first of all evident that if , 



R 



2 » '** , R n are known, then Rj , R^ , ••• , are uniquely determined 



by the p. . distributions. Where p, . is the distribution of material 
y 1J 
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from the i^ initial stage to the final stage. Since it must hold 

that 



j = 1 



p.. =1, i - 1,2, 



and 



n 



R I = p il R i 

1=1 



R II P i2 R i 

1 = 1 




n 



i = 1 



P iP R i* 




Also, it is not difficult to see that the method of solution presented 
for n = 2 earlier becomes increasingly complex as n increases. The 



27 



true computed mean output rate that is obtained by applying the 
approximate procedure, is found, as before, by determining the 
combination of R^+R 2 + +R n = R that yields maximum R. However, 
it would seem that an iterative procedure could be developed to solve 
the problem efficiently for n = 3 or n = 4 on todays high speed 
computers. Efficient extension beyond this point would probably 
require knowledge the R^'s explicitly in terms of the service rates, 
queue lengths and flow distributions. 
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APPENDIX A 




(Distribution from (3,1) to (*4,1)) 
Values of other parameters : 

^ 11 »^ 12 »^ 21»^31 = 0,5 ^ 22*^32 = 0,14 
^41*^42*^13 = 1/3 M 21, M 22, M 31, M 4l, M 42, M 43 = 5 

M 32 =4 qi = 1/3, i = 1,2,3 P 2 = P 3 = (l-Pi)/2 

PARAMETRIC STUDY 1. 



29 






APPENDIX A 



M 

e 

a 

n 

0 

u 

t 

P 

u 

t 

R 

a 

t 

e 

1 

0 
a 
d 
s 
/ 
m 

1 
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[ - Rn and R m 




.3- 





.3 



-h 



-f* H 

.5 .6 



! 1 h 

.7 .8 .9 1.0 



(Distribution from (3,1) to (^,1)) 



Values other parameters : 

A2\'M 31 = 0,f A 22 */^32 c ’ 

A\\\ t A\\2*MkZ ~ 1/3 M21» M 22» M 31> Ml 41» Ml J2» Ml J3 = 5 

M^ 2 « P 2 = P 3 = (l-p^/2 q 1 = 0.2 q 2 « » 0.H 



PARAMETRIC STUDY 2. 
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APPENDIX A 




p-^ (Distribution from (3,1) to (4,1)) 

Values of other parameters ; 

^11 »/^12 *^ 21*^31 = /^ 22V^32 = 

42* / ^43 = ^ 21 *^ 22 *^ 31 *^ 41 *^ 42*^43 = ^ 

M 32 = H P 2 = P 3 = (1 ” p i )/2 q i = 0414 q 2 = q 3 = °* 3 

PARAMETRIC STUDY 3. 
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Values of other parameters : 

Pl “ Qi - 1/3, 1 - 1,2,3 



y^l4l»y^l42»/%3 = 
5 



PARAMETRIC STUDY 
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Values of other parameters : 

21^22V^31^32 = 0,5 

” 1//3 M 21 ,M 22 , %1 ,M 42 , %3 " 5 

Pi = Qi - 1/3, 1 = 1,2,3 



PARAMETRIC STUDY 5. 
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COMPUTER PROGRAM TO OBTAIN APPROXIMATE SOLUTION FOR THE 
MEAN OUTPUT RATE FOR THE REPLENISHMENT AT SEA PROBLEM 

NOTE: TO SIMPLIFY THE COMPUTATIONS THE PROGRAM HAS BEEN 
WRITTEN TO SOLVE ONLY FOR THE CASE WHERE R1 IS GREATER 
THAN OR EQUAL TO R2. THUS, DATA SHOULD BE ENTERED TO 
MAKE THIS SO. INSPECTION OF THE HANDLING RATES FOR REAL 
WORLD DATA WILL NORMALLY INDICATE WHICH PARALLEL PATH 
WILL HAVE THE LARGEST MEAN OUTPUT RATE 

DIMENSIONED VARIABLES 

YY ( I ) VALUES ARE USED AS THE INITIAL SOLUTION FOR THE 
RATIO OF THE EFFECTIVE ARRIVAL RATE TO THE EFFECTIVE 
SERVICE RATE VALUES. VALUES CORRESPOND TO l/IO OF THE 
INTERVAL (0,1) AND DEPEND ON THE APPROXIMATE MEAN 
SERVICE RATE. 

M ( I ) DEFINE QUEUE SPACE AVAILABLE AT FACILITIES (2,1) 
(2,2), (3,1), (3,2), (4,1), (4,2), (4,3). 

X S ( I ) ARE THE PROBABILITIES THAT THE FACILITIES ARE 
BUSY AND IS COMPUTED WITHIN THE PROGRAM FOR FACILITIES 
( 1 , 1 ) - ( 2 , 3 ) . 

U ( I , J ) ARE THE EXPONENTIAL SERVICE TIME INPUTS FOR THE 
FACILITIES. NOTE U ( 1 , 3 ) -U( 3 , 3 ) ARE DUMMIES AND ARE NOT 
USED IN THE COMPUTATIONS. 

P ( I ) IS THE DISTRIBUTION FROM THE UPPER PATH TO THE FINAL 
STAGE . 

Q ( I ) IS THE DISTRIBUTION FROM THE LOWER PATH TO THE FINAL 
STAGE. 

RHO(I) IS AN INTERMEDIATE VALUE USED IN FINAL COMPUTATION. 

RR ( I ) ARE THE MEAN OUTPUT RATES FROM THE FINAL STAGE. 

DIMENSION YY(I0),M(7),XS(2,3),U(4,3),P(3),Q(3) 

DIMENS ION RHO ( 3) ,RR( 3 ) 

WRITE (6,600) 

600 F0RMAT(TI1 ,' WRITE VALUE') 

R E AD ( 5 , 5 ) MM 
5 FORMAT ( 12 ) 

READ( 4,46) YY 

46 FORMAT ( I0F6. 2 ) 

READ(4 ,47 ) M 

47 FORM AT ( 71 6 ) 

RE AD ( 4 ,48 )U 

40 FORMAT ( I2F6. 4 ) 

READ( 4,49)P,Q 
49 FORMAT (6F6 .4 ) 

X S( I , I ) = I . 

X S ( 1 , 2 )= I . 

1 = 1 

DEN= 100. 

AR IS THE INITIAL TRIAL SOLUTION FOR THE MEAN OUTPUT RATE. 

AR=0.6*(U( 4, 1 ) +U( 4, 2 ) + U( 4,3)) 

R I IS THE TRIAL VALUE FOR THE UPPER ROUTE 
Ri=AR/2. 

R 2 IS THE TRIAL VALUE FOR THE LOWER ROUTE 
R 2= AR-R 1 

63 C=Ri/U( 1,1) 
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SUBROUTINE SELEC DETCRMI NES AN INITIAL TRIAL SOLUTION (Y) 
FROM THE TABLE YY FOR THE POLYNOMIAL TO BE SOLVEO BY THE 
NEWTON RAPHSON METHCO. 

CALL SELEC (C ,YY, Y ) 



SUBROUTINE APROX SOLVES THE POLYNOMIAL TO DETERMINE THE 
RATIO OF EFFECTIVE ARRIVAL RATE TO EFFECTIVE SERVICE RATE 
(Y) AND THEN COMPUTES XS(I,J). 

CALL APROX(C,M( 1),Y,XS(2, 1) ) 

C=R2/U ( I » 2 ) 

CALL SELEC (C ,YY,Y) 

CALL APROX(C»M(2) »Y»XS(2»2) ) 

C=R1/(U(2,1)*XS(2,1) > 

CALL SELEC (C ,YY,Y) 

CALL APROX(C»M(3)*Y,XS(3» l) ) 

C=R2/(U(2,2)*XS(2,2) ) 

CALL SELEC(C,YY,Y) 

CALL APROX (C,M<4) , Y,XS{ 3, 2) ) 

A=XS(3,1)*U(3,1) 

B=XS ( 3,2) *U( 3,2) 

DC 60 K= 1 , 3 

R HO ( K ) — ( P C K)*A+Q( K)*B)/U(4,K) 

D=RH0(K)**(M(4+K) + 1) 

60 RR(K)= (RHO(K)-D)/ ( l.-D)*U(4,K ) 

RS IS THE COMPUTED VALUE OF THE OUTPUT RATE. 

RS=RR( 1 ) + RR ( 2 ) +RR I 3 ) 

THIS SECTION OF THE PROGRAM COMPARES THE TRIAL VALUE AR 
TO THE COMPUTED VALUE RS AND ITERATES UNTIL THE DIFFERENCE 
IS LESS THAN .001. THE SECTION ALSO CAUSES ITERATIONS TO 
BE MADE TO FIND THE VALUES OF R1 AND R2 THAT YIELD MAX RS . 

IFII.EQ.I ) GO TO 61 
I FIRS. LT. RSI ) GO TO 62 
RS1=RS 

64 R1=R1+1./DEN 
R 2=AR- R 1 

GO TO 63 

61 R S1 = RS 
1=2 

GO TO 64 

62 R 1=R 1- 1 • / DEN 
R 2 = A R - R 1 

IF! I . NE. 2 ) GO TO 65 
D EN= 1000 . 

1=3 

GO TO 64 

65 1 = 1 
DEN= 100. 

S=R1/R2 

T = ABS I RS-AR) 

I F IT. LE. 0.002) GO TO 66 
IFIRS.GT .AR) GO TO 67 
AP-AR-T/2. 

69 IFIR1.EQ.R2) GO TO 68 
R3=ABS(R1-R2 ) 

IF (R 3-. 1)68,68,70 

70 R 1= AR/ 2 . +R 3- . 1 
R2= AR— R 1 

GC TO 63 

67 AR=AR+T/2. 

GO TO 69 

68 Rl=AR/2. 
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on 



R 2 = AP -R 1 
GC TO 63 

66 WRITE(MM,55)RS»Rl , R2 » RR 

55 FCRMAT(//T5, 'QUEUE OUTPUT I S • ,T2 3 , F 8 . 5// T 5 , • P ARAL LEL 
OUTPUTS ARE' ,T28, 2F8.5//T5, 'TERMINAL OUTPUTS ARE', 

T 2 8 , 3F 8 . 5 ) 

WRITE (MM, 52)U 

52 FORMAT ( //T44 , ' HANOLING RATES •// < 1HO, T 10, 12F7 . 4 ) ) 
WRITEIMM, 53IM 

53 FORMAT ( //T40 » ' MAX I MUM QUEUE L ENGTHS ' //( 1H0, T25, 71 7) ) 
WPITF(MM,54)P,Q 

54 FORMAT (//T37, 'DISTRIBUTION TO FINAL STAGE'// 

( I HO ,T25,6F8.4 ) ) 

S TCP 
END 



S LBROU TINE APROX ( VAL , N, Y , X ) 

AN=N 

30 A = Y * ♦ N 
B = A* Y 

C= ( Y-B )/( l.-B)-VAL*Y 

D=(1.-(AN+1. ) *A+AN*B)/( l.-B)**2-VAL 
Yl=Y-( C/D ) 

IF(ABS(Y1-Y)-. 0001)20,20,22 
22 Y = Y 1 

GO TO 30 

20 X= ( Y-B )/( 1 .-B ) 

RETURN 

END 



SUBROUTINE SELEC(C,A,B) 
DIMENSION A ( 1 0 ) 

DO 25 1=1 , 10 

I F ( C-F LOAT ( I I / 10 • )23 t 23,25 
25 CONTINUE 
2 3 B=A( I ) 

RETURN 

END 



SAMPLE DATA DECK 
C 

YY ( I ) INITIAL TRIAL SOLUTION FOR NEWTON RAPHSON POLY- 
NOMIAL SOLUTION, IT OEPENDS ON THE VALUE OF AR ASSUMED. 

YY ( 1 ) Y Y ( 2 ) YY ( 3 ) YY(4> YY(5) YY(6) YY(7) YY(8) YY ( 9 ) YY(IO) 

20.00 6.00 4.00 2.80 2.20 1.80 1.50 1.20 0.90 0.60 

C 

M( I ) MAX QUEUE LENGTHS 

M2 1 M 22 M3 1 M32 M4 1 M42 M43 

5 5 5 4 5 5 5 



U ( I , 


J) - HANDLING RATES (NOTE 


FUR ACTUAL 


I NPUT 


ALL DATA IS 


ON A 


SINGLE 


CARD WITH 12F6.3 


FORMAT ) . 








U 1 1 


U2 1 


U31 U41 U 1 2 


U22 U31 


U4 2 


U 1 3 


U23 


.5000 


. 5000 


.5000 .3333 .5000 


.4000 .4000 


. 3333 


. 5000 


. 5000 


U33 


U4 3 












. 5333 

r 


. 3333 












p ( i ) 


AND Q( I ) DISTRIBUTIONS 










p( i) 


P( 2) 


P ( 3 ) Q( l) Q ( 2 ) 


Q ( 3 ) 








. 3333 


.3333 


.3333 .3333 .3333 


. 3333 
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